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Abstract 

A wealth of Galactic accreting X-ray binaries have been observed both in low/hard state and high/soft state. The 
transition between these two states was often detected. Observation shows that the transition luminosity between 
these two states is different for different sources, ranging from 1 % to 4% of the Eddington luminosity. Even for 
the same source the transition luminosity at different outbursts is also different. The transition can occur from 
0.0069 to 0.15 Eddington luminosity. To investigate the underlying physics, we study the influence of viscosity 
parameter a on the transition luminosity on the basis of the disk-corona model for black holes. We calculate the 
mass evaporation rate for a wide range of viscosity parameter, 0.1 < a < 0.9. By fitting the numerical results, we 
obtain fitting formulae for both the transition accretion rate and the corresponding radius as a function of a. We find 
that the transition luminosity is very sensitive to the value of a, L/L^dd oc a 2 34 . For 0.1 < a < 0.6, the transition 
luminosity varies by two orders of magnitude, from 0.001 to 0.2 Eddington luminosity. Comparing with observations 
we find that the transition luminosity can be fitted by adjusting the value of a, and the model determined values of 
a are mostly in the range of observationally inferred value. Meanwhile we investigate the truncation of the disk in 
the low/hard state for some luminous sources. Our results are roughly in agreement with the observations. 

Key words: accretion, accretion disks — X-rays: individual (GX 339-4, GS 2000+251) — X-rays: stars 
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1. Introduction 

It is well known that there are two basic spectral states, a 
high/soft state and a low/hard state in accreting X-ray binaries. 
At the high/soft state, the accretion is dominantly via a stan- 
dard disk (Shakura & Sunyaev 1973); while at the low/hard 
state, the accretion is dominantly via an advection dominated 
accretion flow (ADAF) or a radiative inefficient accretion flow 
(RIAF) (Narayan & Yi 1994, 1995a,b; Abramowicz et al. 
1995; for reviews see Narayan 2005 and Kato et al. 2008). 
Similar accretion models are proposed in active galactic nuclei 
(Narayan & Yi 1994, Yuan et al. 2003). Detailed fits to the 
observational spectra show that in the soft/high state of black 
hole the thin disk extends down to the innermost stable cir- 
cular orbit (ISCO); while in the hard/low states, the accretion 
flow consists of two radial zones, i.e., an inner advection dom- 
inated accretion flow (ADAF) that extends from ISCO (or even 
inside the ISCO, see Watarai & Mineshige 2003) to some trun- 
cation radius and an outer thin accretion disk, above which is 
the corona which provides hot mass accretion flow for the inner 
ADAF (Esin et al. 1997). 

The mechanisms to facilitate the state transition have been 
investigated by Honma (1996), Manmoto & Kato (2000) and 
Lu et al. (2004), which are based on a radial conductive en- 
ergy transport process, and by Meyer et al. (2000b), Rozanska 
& Czerny (2000), Spruit & Deufel (2002), and Dullemond & 
Spruit (2005), which are based on a vertical evaporative pro- 
cess. In strong ADAF principle, the transition luminosity is 
defined as the critical accretion rate for an ADAF to exist 



(Narayan &Yi 1995b; Abramowicz et al. 1995) and thus is 
derived from the assumption that all of the accretion energy is 
transferred to electrons and eventually radiates away (<rf s <~ 
<7ie). The disk-corona evaporation model (Meyer & Meyer- 
Hofmeister 1994; Meyer, Liu, Meyer-Hofmeister 2000a; b; 
Liu et al. 2002) can naturally explain both the transition of 
spectral state and truncation of the thin disk in the low/hard 
state. The interaction between the disk and the corona leads to 
mass evaporating from the disk to the corona. The evaporation 
rate reaches a maximum value at a few hundred Schwarzschild 
radii. When the accretion rate feeding the accretion in the outer 
disk is higher than the maximal evaporation rate, the disk ex- 
tends down to the innermost stable orbit. When the accretion 
rate is lower than the maximal evaporation rate, the disk is first 
truncated at the most efficient evaporation region and then the 
inner disk could be completely depleted if the accretion rate is 
too low. The final steady state at high accretion rate is then 
the soft state dominated by disk accretion and at low accretion 
rate the hard state dominated by the ADAF/RIAF. The maximal 
evaporation rate corresponds to the transition accretion rate. 

Investigation of the soft-to-hard state transition for X-ray bi- 
naries with good measurements of flux, distance and mass of 
the compact object shows that the transition luminosity is at 
about 1-4% of the Eddington luminosity, with a mean of 2% 
(Maccarone 2003). This is in good agreement with the theo- 
retical prediction of 2% by the disk corona evaporation model 
(Meyer et al 2000a). Nevertheless, the transition luminosity of 
some sources deviates this value largely. For instance, the tran- 
sition luminosity observed for GX 339-4 reaches 0.15 L^dd, 



2 



Erlin Qiao & B.F. Liu 



[Vol. 



while it is only 0.0069 L E dd for GS 2000+251. Studies of in- 
dividual object, GX 339-4, show that even for the same object 
the transition occurred at different luminosity during different 
outbursts (Zdziarski et al. 2004, Belloni et al. 2006). Even 
in the same outburst, the transition luminosity can differ in the 
rise and decay phases (Miyamoto et al. 1995). If it is the disk 
evaporation that causes the state transition in above cases, there 
should be something changed in the corona. 

One of the potential factors leading to different transition lu- 
minosity is the change of viscosity parameter, which is fixed to 
0.3 in previous studies (Meyer et al. 2000a;b). Although the 
effective temperature profile and thus spectra do not depend 
on the a value at all in the standard disk, the emissivity in an 
optically thin corona depends on the density and the density 
is roughly inversely proportional to a. Therefore, the radia- 
tive luminosity from an optically thin corona/ADAF depends 
on the value of viscosity parameter. On the other hand, there 
is a long research history of the a value through the study of 
outbursts of dwarf novae and X-ray novae (Meyer & Meyer- 
Hofmeister 1983; Smak 1984; Lin et al. 1985; Mineshige & 
Wheeler 1989). The most recent study for black holes (King , 
Pringle & Livio 2007) yields 0.1 < a < 0.6 (if the viscosity is 
defined as v = ^aV s H). 

For the range of 0.1 < a < 0.3, the influence of viscosity 
has been studied by Meyer-Hofmeister & Meyer (2001) for 
a few number of a value. In their calculations, the temper- 
ature of electrons is assumed equal to that of ions, that is a 
one-temperature model. In reality, the accretion flows around 
a black hole can reach a temperature of ~ 10 9 K or higher, at 
this situation electrons and ions decouple. Therefore, we take 
the two-temperature disk corona model (Liu et al. 2002) and 
carry out more detailed calculations for a series of viscosity 
parameters in the range of0.1<a<0.9, which covers the 
recent observational estimates of the viscosity range (King et 
al. 2007). 

In this work, fitting formulae of the numerical results for 
both the maximal mass accretion rate and the corresponding 
radius depending on the viscosity parameter a are given. On 
the basis of this, we compare our theoretical results with ob- 
servations for individual X-ray binary and AGNs. In Sect. 2 we 
briefly describe the physics of the corona above the thin disk 
and list equations. In Sect.3 we show our numerical results and 
in Sect.4 we compare the model predictions with observations. 
Our discussions and conclusion are given in Sect.5 and Sect. 6. 

2. The model 

The disk corona model for accreting black holes is estab- 
lished in detail by Meyer et al. (2000a) and later this model is 
extended close to the central black hole by taking into account 
the decoupling of ions and electrons (Liu et al. 2002). It is 
further developed by including the inflow and outflow of mass, 
energy and angular momentum from and towards neighboring 
zones (Meyer-Hofmeister & Meyer 2003). The basic physics 
of the model is briefly described below. 

We consider a hot corona above a geometrically thin stan- 
dard disk around a central black hole. In the corona, viscous 
dissipation leads to ion heating, which is partially transferred 
to the electrons by means of Coulomb collisions. This energy 



is then conducted down into lower, cooler and denser corona. 
If the density in this layer is sufficiently high, the conductive 
flux is radiated away. If the density is too low to efficiently 
radiate the energy, cool matter is heated up and evaporation 
into the corona takes place. The mass evaporation goes on 
until an equilibrium density is established. The gas evaporat- 
ing into the corona still retains angular momentum and will 
differentially rotate around the central object. By friction the 
gas looses angular momentum and drifts inward thus contin- 
uously drains mass from the corona towards the central ob- 
ject. This is compensated by a steady mass evaporation flow 
from the underlying disk. The process is driven by the gravita- 
tional potential energy released by friction in the form of heat 
in the corona. Therefore, mass is accreted to the central ob- 
ject partially through the corona (evaporated part) and partially 
through the disk (the left part of the supplying mass). In the 
inner region the evaporation becomes so efficient that all the 
matter accreting through the disk is heated up and accreted to- 
wards the black hole via the corona. Only when the accretion 
rate is higher than the evaporation rate, can the disk survive. 
Therefore, the disk is truncated at the radius where the evapo- 
ration rate equals to the accretion rate. In the case of accretion 
rate higher than the maximal evaporation rate, the disk extends 
down to the ISCO. In this work, we investigate how the maxi- 
mal evaporation rate and the corresponding radius varies with 
viscosity parameter a. We use the model of Liu et al. (2002) 
and include the modification on the energy equation. For clar- 
ity, we list the equations describing the physics of corona as 
follows: 

Equation of state 



(i) 



where /i = 0.62 is the molecular weight assuming a standard 
chemical composition (X = 0.75, Y = 0.25) for the corona. For 
convenience, we assume the number density of ion m equals 
to that of electron n e , which is strictly true only for a pure 
hydrogen plasma. 

Equation of continuity 
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Equation of the z-component of momentum 
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where tjm is the mass advection modification term and t\e is 
the energy modification term. We take tjm = 1 for the case 
without consideration of the effect of mass inflow and out- 
flow from and towads neighboring zones in the corona, and 
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Ve = f]u + 0.5 is a modification to previous energy equations 
(for details see Meyer-Hofmeister & Meyer 2003). qi e is the 
exchange rate of energy between electrons and ions, 
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where m p and m e is the proton and electron mass, k is 
Boltmann constant, a is the Thomas scattering cross and In A = 
20 is the Coulomb logarithm. 

The energy equation for both the ions and electrons 
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where n e n,i(T) is the bremsstrahlung cooling rate and F c 
the thermal conduction (Spitzer 1962), 



dz 



(8) 



with «o = 10 _6 ergs _1 cm _1 K -7 / 2 for fully ionized plasma. 
All other parameters in above equations are under standard def- 
inition, and are in cgs units. The five differential equations, 
Eqs.(2), (3), (4), (7) and (8), which contain five variables P(z), 
T e , F c , and m(z)(= pv z ), can be solved with five boundary 
conditions. 

At the lower boundary zq (the interface of disk and corona), 
the temperature of the gas should be the effective temperature 
of the accretion disk. Previous investigations (Liu, Meyer, & 
Meyer-Hofmeister 1995) show that the coronal temperature in- 
creases from effective temperature to 10 6 5 K in a very thin 
layer and thus the lower boundary conditions can be reason- 
ably approximated (Meyer et al. 2000a) as, 

Ti = T e = 10 6 5 K, and F c = -2.73 x 10 6 P at z = z . (9) 

At infinity, there is no pressure and no heat flux. This requires 
sound transition at some height z=z\. We then constrain the 
upper boundary as, 

sp 

F c = and v 2 = V 2 = P/p = —{T l + T e ) at z = z v (10) 

With such boundary conditions, we assume a set of lower 
boundary values for P and m to start the integration along z. 
Only when the trial values for P and m fulfill the upper bound- 
ary conditions, can the presumed P and m be taken as true 
solutions of the differential equations. 

Here we wish to point out that the viscosity parameter a we 
are studying in this work only appears in the viscous heating 
rate, ^aPQ, and radial drift speed, vr oc a, which is eventually 
associated with the heating and cooling rates. 



3. Numerical results 

In our calculation we fix the mass of black hole as M = 6M Q 
and vary the viscous parameter a in the range of 0.1-0.9. For 
every given a, the evaporation rate is calculated as m = pv z at 
the lower boundary and then is integrated in the radial one-zone 
region, M cvap w 2-kR 2 vtlq. This integrated value represents the 
evaporation rate in a given region around distance R and thus 
can be compared with the accretion rate. 

Fig. 1 shows how the evaporation rate varies with the distance 
from the black hole. Here the evaporation rate is scaled by the 
Eddingtonrate, M Edd = L Edd /rjc 2 = 1.39 x 10 18 M/M Q gs _1 
(where ?/ = 0.1 is the energy conversion efficiency) and the ra- 
dius is scaled by the Schwarzschild radius, Rg = 2GM/c 2 = 
2.95 x 10 5 M/Mq cm. It can be seen that for every given a 
the evaporation rate reaches a maximum at a certain distance 
and decreases very quickly outside the maximal evaporation 
region. The occurrence of maximal evaporation rate is a result 
of energy balance in the disk corona system. In the outer re- 
gion, the evaporation rate increases toward central black hole 
since the released accretion energy increases (oc i? _1 ), which 
is the source of energy for evaporation. However, with increase 
of number density in the corona, the radiation becomes more 
and more efficient (since the bremsstrahlung is proportional to 
the square of number density) and only a very little (or even 
no) gas is heated up and evaporates. Thus, the evaporation rate 
reaches a maximum value and decreases again in the inner re- 
gion. 

Fig. 1 also shows that the value of evaporation rates increases 
with the value of viscosity parameter, which is significant in the 
inner region but not obvious in the outer region. The increase 
of evaporation rate with viscosity is understood as a result of 
efficient heating and unimportant advection. In the inner re- 
gion, the advective cooling (related to a) is not the dominant 
cooling process, as shown by our numerical data. For larger 
value of a, the viscous heating is more efficient, which leads 
to more energy that is conducted down to the chromosphere 
above the disk. There the conducted energy is partially radi- 
ated away and partially used to evaporate the cool gas. Thus, 
the evaporation rate enhances for larger a. This effect is not 
significant in regions far away from the black hole since most 
of the viscous heating is transferred to the internal energy of 
corona gas, only a very small part is conducted down and rele- 
vant to evaporation. 

The maximal evaporation rate is very sensitive to the value 
of a. The values of the maximal evaporation rate for a=0.1 
and a=0.9 are 0.001 and 0.2 Eddington rate, respectively, in- 
creasing by a factor of about 200 from a = 0.1toa = 0.9. The 
location of the maximal evaporation rates for a=Q. 1 and a=0.9 
are ~ 1880i?s an d 25i?s, decreasing by a factor of about 75 
from a = 0.1 to a = 0.9. The detailed results associated with 
the maximal evaporation rate are listed Table 1 . 

Plotting the data in Fig. 2, we find that the dependence of 
both the maximal evaporation rate and the corresponding ra- 
dius on the viscous parameter can be expressed in power-law 
forms. The best fits to the data are given as 



(M evap /M E dd) r 



0.38a 2 ' 34 , 



(11) 
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Fig. 1. The distribution of evaporation rate along distances for a series 
of viscosity parameters. From bottom to top a is 0.1, 0.15, 0.2, 0.3, 
0.5, 0.9 respectively. It shows that the evaporation rate increases with 
increasing a. For larger a the maximal evaporation rate is higher and 
the corresponding radius is smaller. 



Table 1. The evaporation feature for different value of a 





a 


^max / Rs 


M max / Mgdd 


M = 6M Q 


0.10 


1883.65 


1.19 x Kr a 




0.15 


891.25 


4.58 x 10~ 3 




0.20 


446.68 


1.01 x 10~ 2 




0.25 


281.84 


1.78 x 10~ 2 




0.30 


223.87 


2.73 x 10^2 




0.40 


112.20 


5.16 x 10~ 2 




0.50 


70.79 


8.11 x 10" 2 




0.70 


37.58 


1.49 x lO" 1 




0.90 


25.12 


2.24 x lO" 1 


M = 10 B M Q 


0.3 


223.87 


2.73 x lO - * 




0.5 


70.79 


8.11 x 10" 2 



Note:M max and R max represent the maximal evaporation rate 
and the corresponding radius for different a. 

(R/Rs)m=m^ » 18.80a" 2 00 . (12) 

As it has been pointed out before (Meyer et al. 2000a;b), the 
existence of a maximal evaporation rate along radial direction 
has very important consequence on the spectral state transition. 
When the accretion rate supplying to the outer disk is lower 
than this value, the disk is completely depleted first at the most 
efficient evaporation region, i.e. the region where the maximal 
evaporation rate reaches. The disk is thus truncated and then 
recedes outwards until the evaporation is too weak to compete 
with the mass flowing rate in the disk. At the final steady state 
the disk is truncated at a distance where the supplying accretion 
rate is equal to the evaporation rate. In the inner region a pure 
corona/ADAF fills in. Therefore, the accretion is via an inner 
ADAF and an outer thin disk, and the radiation from the ac- 
cretion flows forms a hard-state spectrum. When the accretion 
rate is higher than the maximal evaporation rate, evaporation 
cannot deplete the disk at any distances, and the disk extends 



down to the innermost stable orbit. Therefore, the accretion is 
dominated by the thin disk, and the corona is very weak. The 
spectrum is dominated by multi-color blackbody. Therefore, 
the value of maximal evaporation rate represents the critical 
mass accretion rate for the transition from hard/low state to 
high/soft state. The corresponding radius indicates how close 
to the black hole the thin disk is before the state transition. 

Previous study (Meyer & Meyer-Hofmeister 2000a, Liu et 
al. 2002) reveals that the properties of evaporation are inde- 
pendent on the black hole mass. To clarify this for the case 
of different viscous parameters, we also perform calculations 
with a black hole mass of 10 8 M Q . We find that for any given 
a in the range of a = 0.1 to 0.9, the results are mass indepen- 
dent as the evaporation rate is scaled in Eddington accretion 
rate and the distance in Schwarzschild radius. As an example, 
we list the value of a = 0.3,0.5, one can see from Table 1. This 
indicates that our model and predictions are not only valid in 
stellar mass black hole, but also are applicable to supermassive 
black holes like AGNs. The most direct implication is that the 
disk truncation and state transition should occur at the same 
accretion rate for both the black hole X-ray binaries and AGNs 
if the values of a are the same. 

4. Comparison with observation 

Our calculations show that for standard viscosity, a = 0.3, 
the maximal evaporation rate is of the Eddington rate, 

which implies the transition from hard to soft state occurs at 
an accretion rate ~2% of the Eddington rate. This is in good 
agreement with Maccarone (2003)'s observational study on the 
X-ray binaries, where it is found the average transition rate is 
also 2% of the Eddington rate. 

The change of viscosity provides a possible explanation for 
the deviation of the transition accretion rate from the average 
value. For X-ray binaries, King et al. (2007) found a range 
from 0.2 to 0.4, which is 0.3 to 0.6 in our definition of a. In 
this range the disk evaporation model predicts a range for the 
transition accretion rate of 2% to 11%. The outer disk can 
reach down to 52i?s before transitions to the soft state. For 
a lower viscous parameter the transition accretion rate can be 
even smaller. For instance, in the case of a = 0.2, the pre- 
dicted accretion rate is 1%. The strong dependence of maximal 
evaporation rate on the viscous parameter indicates that a small 
fluctuation of viscosity could lead to the observational inferred 
variation in the transition rate and may explain why the spectral 
state transits from hard to soft state at different accretion rates 
for different objects, or for the same object at different out- 
bursts. The dependence of the truncation radius at state transi- 
tion on the viscous parameters also provides a possible mech- 
anism for the change of observational inferred disk truncation 
extending down to ~ 100i?s at hard state before transition . 
The numerical fitting formulae given in this work are useful 
to determine the value of the viscosity parameter from the ob- 
served transition luminosity and then determine the truncating 
radius for individual objects. On the other hand, if viscosity 
can be determined by some other way, Eqs.(ll) and (12) can 
predict the accretion rate at spectral transition and spectral en- 
ergy distribution contributed by the disk and the corona. 
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Fig. 2. LeftThe maximal evaporation rate vs. viscosity parameter. Right:The radius corresponding maximal evaporation rate vs. viscosity parameter. 
The figure shows that the maximal evaporation rate increases and the corresponding radius decreases with increase of viscous parameter a 



4.1. Fits to the transition luminosity 

From observations one obtains the luminosity at the spectral 
transition. Adopting r\ = 0.1 for the energy conversion coef- 
ficient, we have itran/^Edd = M I M-Edd- Thus, according to 
our evaporation model the transition luminosity is a function of 
viscous parameter a, that is, 

£tra„/£ E dd« 0.38a 2 34 . (13) 

To fit the transition luminosity, we collect data from literatures 
for the sources that the transition luminosity is well measured. 
The detailed results are listed in Table 2. If the variation of 
transition luminosity is caused by the change of viscosity, we 
can calculate the value of a from equation (13). The derived 
values are also listed in Table 2. One can see that though the 
transition luminosity varies largely from about 0.0069 L^dd in 
GS 2000+251 to 0.15 L E dd in GX 339-4, the values of a are 
in a relatively small range from 0.18 to 0.67. Here we must 
note the different definition of a in different literatures. We 
take v = ^aVgCl (Shakura & Sunyaev 1973) for the kinematic 
viscosity, while King et al. (2007) adopt v = aV/O (Frank 
et al. 2002). Thus, the value of a in our work differs by a 
factor of 3/2 from that of King et al. (2007). Translating the 
above values of viscous parameters to the definition of King 
et al. (2007), the fitted value of o. is in the range of o. ~ 0.12 
- 0.45. These values are in agreement with the observational 
values of 0.1- 0.4 inferred by King et al. (2007). Here we 
wish to point out that there are uncertainties in the inferred 
transition luminosity Ltran/^Edd due to the uncertainties in 
measuring the mass of black hole and its distance. This leads 
to a small deviation in the value of viscosity parameter in fitting 
the transition luminosity. Nevertheless, we expect the value of 
a is still in the range of King et al. (2007). 

With the value of viscous parameter a, we can calculate 
the truncation radius of the disk at the spectral transition from 
equation (12). If observations can determine the transition ra- 
dius, say, by means of emission lines, one can compare the ob- 
servationally inferred value with the model prediction, thereby 



being able to check whether the transition luminosity is caused 
by the change of the viscosity. 

4.2. Fits to the truncation radius in the low/hard state 

In the low/hard state of black hole accreting systems, the ac- 
cretion rate is smaller than the maximal evaporation rate, thus 
the disk is truncated at a certain distance where the accretion 
rate is equal to the evaporation rate. For a standard viscos- 
ity parameter, a = 0.3, the evaporation model predicts that the 
truncation radius is larger than 224i?s, as listed in Table 1. This 
value is larger than the values inferred from observations for 
some objects. Previous investigation on the effect of magnetic 
fields and conduction (Qian et al. 2007; Meyer-Hofmeister & 
Meyer 2006) show that the truncation can shift to ~ 100i?s if 
the magnetic field and reduced conduction are taken into ac- 
count. However, for some objects, this truncation radius is still 
too large to account for the observation. This contradiction can 
be alleviated by taking into account the fluctuation of the vis- 
cous parameter a. Since the truncation radius in systems with 
accretion rate less than ~ 2-3 percent of Eddington rate can be 
well explained by changing the strength of magnetic field and 
thermal conduction for standard viscous parameter (Qian et al. 
2007), here we only collect observational data with the mass 
accretion rate larger than this value. 

NGC 4593 is known as a Seyfert galaxy. By fitting the con- 
tinuum spectrum of optical-UV and Fe K a line profile, Lu & 
Wang (2000) find the observational evidence supporting that 
thin disk truncates at > 30i? s . The corresponding mass ac- 
cretion rate in the outer thin disk is about 0.0557\/Edd- This 
requires that a value of a > 0.44 from equation (11), with 
which the maximal truncation radius is about 98 i?s according 
to equation (12). If a large a is associated with strong mag- 
netic fields, the strong magnetic fields in this object will result 
in a truncation radius smaller than 98 i?s- The discrepancy of 
a factor of 3 between the observation and theoretical predic- 
tion could be decreased by taking both the magnetic field and 
thermal conduction into account (Qian et al. 2007). Therefore, 
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we expect that the very small truncation of the thin disk in the 
low/hard state of NGC 4593 is a combined effect of viscos- 
ity, magnetic field and conduction. Note that the viscosity and 
conduction are mutually related to the magnetic fields. 

GX 339-4 is of well known black hole X-ray binaries. 
Zdziarski et al. (2004) study the 15 outbursts of GX 339-4 from 
1987 to 2004. When in the low/hard state the highest mass ac- 
cretion rate can reach — 0.26 MEdd> an d also at different time 
the accretion rate can be — 0.14, 0.17, 0.07, 0.02 M E dd in the 
low/hard state (Zdziarski et al. (2004). From equation (1 1) this 
requires the value of a larger than 0.85, 0.65, 0.71, 0.49, 0.28 
respectively. For these value of a, the maximal truncation ra- 
dius is derived from equation (12) as 26, 44, 37, 80, 231 i?s- 
These results are in agreement with the observational values of 
Ri n —20-200 Rs, although the discrepancy remains. 

5. Discussion 

The disk-corona evaporation model provides a physical 
mechanism for the spectral state transition between low/hard 
state and high/soft state and for truncation of the thin disk in the 
low state. Nevertheless, it should be pointed out that there are 
still uncertainties in quantitative fits to observations. The value 
of the mass evaporation rate in our model is affected by several 
factors, such as the magnetic field (Qian et al. 2007), the ther- 
mal conduction (Meyer-Hofmeister & Meyer 2006), and the 
inverse Compton scattering (Liu et al. 2002). It also depends 
sensitively on the viscosity parameter a, as shown from our 
calculations. Therefore, one should take into account all these 
factors in quantitative fits to observations. 

5.1. The effect of Compton cooling 

Compton upscattering of photons from the disk can be im- 
portant in cooling of the corona (Liu et al. (2002). At high ac- 
cretion rates, the thin disk extends down to the ISCO, Compton 
cooling becomes efficient and hence less heat is conducted 
downwards, which leads to a decreased evaporation rate. The 
maximal evaporation rate also decreases, indicating that the 
transition luminosity from soft to hard decreases as a result of 
Compton cooling. On the other hand, at low accretion rate the 
disk is truncated. Compton upscattering of disk radiation is not 
efficient, thereby hardly affects the evaporation rate. Therefore, 
our model in this work is valid for the transition from low/hard 
state to high/state. For the transition from high/soft to low/hard 
state, the transition luminosity is a few times lower than that 
from low/hard to high/soft transition (Meyer-Hofmeister et al. 
2005; Liu et al. 2005), producing the observed hysteresis. It 
is interesting to note that different alpha values may also con- 
tribute to the hysteresis. This makes our prediction on hystere- 
sis closer to the observations for objects showing very large 
hysteresis (e.g. Miyamoto et al 1995). 

In the case of a weak inner disk existing at low/hard states 
(Liu et al. 2007), the Compton effect can also lead to a lower 
evaporation rate, but it can not change the maximal evaporation 
rate significantly. 



5.2. The combined effects of magnetic fields, conduction and 
viscosity 

As it has been shown in previous study (Qian et al. 2007), 
the magnetic fields and conduction hardly affect the transition 
luminosity, but lead to a small truncation radius at the state 
transition. Here it is shown that increase of viscosity results 
in not only a decreased transition radius but also an enhanced 
transition luminosity. Combination of these effects can make 
our model more flexible in fits to the observationally inferred 
transition luminosity and and the truncation radius. Here one 
should note the physical link between magnetic field strengths 
and viscosity parameter. In fact, the viscosity is mostly pro- 
vided by magnetic stress, though we are not very clear the de- 
tailed dependence of a on magnetic stress. 

5.3. The value of viscous parameter a 

Collecting and analyzing observational results from dwarf 
novae, soft X-ray transients and AGNs, King et al.(2007) find 
the range of viscous parameter is 0.1 < a < 0.4. Fits to the 
observational transition luminosity by our model give similar 
range for the viscous parameter. However, the value of a de- 
rived by many MHD simulations is less than 0.02. For instance, 
Stone et al. (1996) derive a <0.01; Brandenburg et al.(1995) 
get 0.001 < a <0.005, Hirose, Krolik & Stone (2006) obtain 
a ~ 0.016, Hawley &Krolik (2001) get a - 0.1. The analysis 
of detailed observational data (King et al. 2007) indicates that 
most MHD simulations underestimate the value of a. 

Caution should be taken to the fact that the a value obtained 
from the numerical simulations depends how to calculate them, 
(e.g., Machida, Hayashi & Matsumoto 2000, Machida & 
Matsumoto 2003, Machida, Nakamura & Matsumoto 2004). 
This is because the magnetized accretion flow shows highly in- 
homogeneous structure and significant time fluctuations. Thus, 
the derived alpha values depend how to make time and spatial 
averages of physical quantities. Generally speaking, the ra- 
tio of magnetic energy to gas energy and thus the alpha value 
increases with increase of height; the alpha values sometimes 
even exceed unity in corona (see, e.g., Machida et al. 2000). 
Hence, if we make the spatial average with weight of magnetic 
energy instead of gas density or gas energy, the derived alpha 
values tend to be higher. 

5.4. The dependence of wind loss on a 

The fraction of mass carried away by the wind sensitively 
depends on the value of a. For a = 0.3, the fraction increases 
from —7% at the maximal evaporation region (R = 10 2 35 i?s) 
to —9.8% at R = 10 4 25 i?s- For a = 0.1, the fraction increases 
from —1.4% at the maximal evaporation region (R~ 10 3 27 i?s) 
to —2.7% at R = 10 4 25 i?s- For a = 0.9, the fraction increases 
from -16% at R = 25R S to 19% at R = 10 4 25 R S . One can 
see the fraction of mass carried away by the wind increases 
systemically with an increase of the value of a. This can be 
understood as the higher temperature and pressure in the larger 
a case. 

5.5. A mechanism for the "strong ADAF" principle 

Our model indicates that the accretion flow changes from a 
corona/ ADAF to a thin disk when the accretion rate exceeds 
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Table 2. Transition luminosity and fitting value of the viscosity parameter in accreting X-ray bina- 
ries 



source 


Aran/ -^Edd 


a 


references 


XTE J 1859+226 


0.12 




0.61 


(Gierlinski & Newton 2006) 


XTE J 1739-278 


0.13 




0.63 


(Gierlinski & Newton 2006) 


XTE J 1650-500 


0.02 




0.28 


(Gierlinski & Newton 2006) 


4U 1543-47 


0.07 




0.49 


(Gierlinski & Newton 2006) 


Nova Mus91 


0.031* 




0.34 


(Maccarone 2003) 


GS 2000+251 


0.0069* 




0.18 


(Maccarone 2003) 


CygX-1 


0.028* 




0.33 


(Maccarone 2003) 


GRO J 1655-40 


0.0095* 




0.21 


(Maccarone 2003) 


LMC X-3 


0.014* 




0.24 


(Maccarone 2003) 


Aql X-l 


0.019* 




0.28 


(Maccarone 2003) 


4U 1608-52 


0.042* 




0.39 


(Maccarone 2003) 


4U 1728-34 


0.050* 




0.42 


(Maccarone 2003) 


XTE J 1550-564 


0.13 




0.63 


(Gierlinski & Newton 2006) 


XTE J 1550-564 


0.03 




0.34 


(Gierlifiski& Newton 2006) 


XTE J 1550-564 


0.034* 




0.36 


(Maccarone 2003) 


GX 339-4 


0.15 




0.67 


(Gierlinski & Newton 2006) 


GX 339-4 


0.06 




0.45 


(Gierlinski & Newton 2006) 



Note: (1) Aql X-l, 4U 1608-52, 4U 1728-34 are neutron stars, the others are 
block holes or black hole candidates. (2) The luminosity band in Gierlinski & 
Newton (2006) is 1.5-12 kev. (3) Asterisk denotes the transition luminosity from 
high/soft state to low/hard state, for which the fitting value of a should be slightly 
higher if the Compton cooling by disk photons is included in our model. (4) 
Different transition luminosities for XTE J 1550-564 and GX 339-4 are from different 
outbursts. 



the maximal evaporation rate. Thus, the maximal evaporation 
rate represents the critical accretion rate for the accretion tak- 
ing the form of a corona/ADAF. Our detailed calculations show 
that the maximal evaporation rate depends on the viscous pa- 
rameter approximately in the form of to cx a 2 34 . This result is 
similar to the critical accretion rate for an ADAF to exist, that 
is, m cr i t oc a 2 (Narayan &Yi 1995b; Abramowicz et al. 1995; 
Mahadevan 1997). The similar dependent laws on the viscos- 
ity are easy to understand since both the maximal evaporation 
rate and critical accretion rate for an ADAF to exist are derived 
from the assumption that all of the viscous heating is trans- 
ferred to electrons and eventually radiated away (q^ is ~ q lc , for 
details see Mahadevan 1997). Now that the evaporation as a 
consequence of disk corona interaction can completely deplete 
the inner disk and turn it into a pure corona/ADAF when the 
accretion rate is less than the maximal evaporation rate, the 
disk evaporation model provides a naturally physical mecha- 
nism for the "strong ADAF" principle, that is, whenever the 
accretion rate is lower than the critical accretion rate, the ac- 
cretion is via an ADAF. 

6. Conclusion 

Using the two-temperature disk corona model, we inves- 
tigate the influence of viscosity on the evaporation feature. 
We find that the evaporation rate depends strongly on the 
value of viscous parameter. For 0.1 < a < 0.6, the maximal 
evaporation rate varies by two orders of magnitude, from 
0.001 to 0.2 Eddington rate. The distance corresponding to 
the maximal evaporation rate also varies by near two orders of 



magnitude. An numerical fitting formulae of the maximal mass 
evaporation rate as a function of viscosity is given by fitting 
the numerical data. The corresponding radius is also given as 
a function of a. Since the maximal evaporation rate represents 
the critical accretion rate at the transition from low/hard to 
high/soft state, we compare our model with observations. We 
find that the transition luminosity can be fitted by adjusting 
the value of a, and the model determined values of a are 
mostly in the range of observationally inferred value (King et 
al. 2007). Meanwhile we investigate the truncation of the disk 
in the low/hard state for some luminous sources. Our results 
are roughly in agreement with the observations. 

We thank F. Meyer and E. Meyer-Hofmeister for their 
comments. This work is supported by the the National 
Natural Science Foundation of China (grants 10533050 and 
10773028). 
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